Neutron star interiors and topology change 
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Quark matter is believed to exist in the center of neutron stars. A combined model consisting of 
quark matter and ordinary matter is used to show that the extreme conditions existing in the center 
could result in a topology change, that is, in the formation of wormholes. 
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Wormholes are handles or tunnels in the spacetime topol- 
ogy connecting different parts of our Universe or of dif- 
ferent universes. The meticulous analysis in Ref. [l[ has 
shown that such wormholes may be actual physical ob- 
jects permitting two-way passage. 

Also introduced in Ref. pj is the metric for a static 
spherically symmetric wormhole, namely 

d s 2 = -e^dt 2 + e x ^dr 2 + r 2 {d0 2 + sin 2 dip 2 ), (1) 

where e A(V ) = 1/(1 — b(r)/r). Here b = b(r) is the shape 
function and (f> — <f)(r) is the redshift function, which 
must be everywhere finite to prevent an event horizon. 
For the shape function we must have b(r ) — r , where 
r = r is the radius of the throat of the wormhole; also, 
b'( r o) < 1 and b(r) < r to satisfy the flare-out condition 
[l[. We are using geometrized units where G = c = 1. 

While Morris and Thorne concentrated on traversable 
wormholes suitable for humanoid travelers (and possible 
construction by an advanced civilization) , a more general 
question is the possible existence of naturally occurring 
wormholes 2]. It has also been suggested that the ex- 
treme conditions inside a massive compact star may re- 
sult in a topology change 3] . It is proposed in this paper 
that a neutron star is a candidate for such a topology 
change as a result of quark matter [3-[6j], which is be- 
lieved to exist at the center of neutron stars 0] ■ To that 
end, we will use a two-fluid model consisting of ordinary 
matter and quark matter based on the MIT bag model 
Q. Two cases will be considered: the non-interacting 
and interacting two-fluid models. (We will concentrate 
mainly on the former case.) The reason for having two 
cases is that quark matter is a viable candidate for dark 
matter and so may be weakly interacting [9(. 

A class of exact solutions of the Einstein-Maxwell 
field equations describing a wormhole with an anisotropic 
matter distribution is discussed in Ref. 

A combined model consisting of neutron-star matter 
and of a phantom/ghost scalar field yielding a nontrivial 
topology is discussed in Ref. [ll[. The model is referred 
to as a neutron-star-plus- wormhole system. 



For our basic equations we follow Ref. [12j. The energy 
momentum tensor of the two-fluid model is given by 

Tq = /^effective = P + Pq, (2) 
T l = T 2 = -^effective = ~(P + Pq)- (3) 

In Eqs. ([2]) and ([3]), p and p correspond to the respec- 
tive energy density and pressure of the baryonic matter, 
while p q and p q correspond to the energy density and 
pressure of the quark matter. The left-hand sides are the 
effective energy density and pressure, respectively, of the 
composition. 

The Einstein's field equations are listed next [l2j |: 
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A' 1 \ 1 „ , 

-) + — = ^( P + p q ), (4) 

r r A J r z 



i) -4 = 87r(p + p g ), (5) 
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iTT(p + Pq ), (6) 



recalling that T/ = T 2 2 . 

In the MIT bag model, the matter equation of state is 
given by 



P q = 7r(Pg-4B), 



(7) 



where B is the bag constant, which we take to be 
145 MeV/ (fm) 3 (slli^. Since quarks are part of the stan- 
dard model, we assume that p q > 0; hence Pg> 4-B. For 
normal matter we use the equation of state [10( 



p = mp, < m < 1. 



(8) 



Since we are assuming the pressure to be isotropic, the 
conservation equation is [121 ] 
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j r -<p Ineffective + ^effective] 

dr 2 



0. 



(9) 
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III. SOLUTIONS OF BASIC EQUATIONS 

As noted in the Introduction, we will concentrate mainly 
on the non-interacting fluid model. This means that the 
two fluids, normal matter and quark matter, do not inter- 
act. The resulting conservation equations are therefore 
independent of each other. Using Eqs. (|7J| and ©, we 
have 



and 



dp _j_ ,/ / 1 + m\ _ ( , 
dr \ 2m 



^ + 24>'( Pq -B)=Q. 



The solutions are 



(l+m)/2m 



and 



Pq = B + p {qfi) e 



-2(4 



(10) 



(11) 



(12) 



(13) 



where po and P( g ,o) are integration constants. Equations 
dU), ©, and © now give 

2 2 r 

87r[p(l + 3m) + (2p g -4S)] (14) 

It becomes apparent that this equation is linear in e~ x 
once it is rewritten in the following form: 



= -(87r)Kl + 3m) + (2p,-4B)]. 
The integrating factor is 

I F = e *+21n0'+41nr 

leading at once to 
d 



dr 



eV)V (Sn)[p(l + 3m) + (2 Pq AB)]. 

To write the solution for e~ A , we first define F{r) = 
p(l + 3m) + (2p q — AB), which becomes, after substituting 
the solutions p and p q , 

F(r) = Po e-^ 1+m ^ 2m {\ + 3m) + 2p {qfi) e-^ - 2B. 

We now have 

16tt 



6') 2 r 4 



e' t '(j)'r\F{r l )dr 1 . (15) 



The lower limit r = ro is the radius of the throat of the 
wormhole, as we will see. 

As noted in Sec. HI the shape function b = b(r) is 
obtained from e~ x ( r \ so that 



(r) = r(l - e" A «) 



We obtain 
b(r) = i 



16vr 



e^0V 4 F(ri)dr 1 



(16) 



(17) 



Note especially that b(r ) — r , which characterizes the 
throat of the wormhole. 

To study the other requirement, b'(ro) < 1, we start 
with 



b'(r) = 1 - e - A(r) + 
which leads to 



,-A(r) 



dr 



b'(r ) = l + r - 



16tt 



-poe _,/,(ro)(1+m)/2m (l + 3m) 
-2 P(q ^e- 2 ^ +2b\ . (18) 



Given that (f>' , the gradient of the redshift function, is 
related to the tidal force, one would expect its magnitude 
to be extremely large. So b'(ro) will indeed be less than 
unity if the expression inside the brackets in Eq. (ITBl is 
negative and sufficiently small in absolute value. If this 
is the case, then the flare-out condition is met, that is, 



r b'(r ) - b(r ) 
2[b(roW 



< 0. 



(19) 



Since it indicates a violation of the weak energy con- 
dition, the flare-out condition is theprimary prerequi- 
site for the existence of wormholes [1|. (The condition 
b(r)/r < 1 will be checked later.) 

Unfortunately, the constants B and P( q .o) are also 
large, so that Eq, (fT8|) must be examined more closely. 



IV. THE GRADIENT OF THE REDSHIFT 
FUNCTION 

To analyze Eq. (fT5|) , we start with the gradient of <p [li 



dr 



2Gm{r) 



4 G ' 
inr—p 



1 - 2Gm(r)/cV 



where m(r) is the total mass-energy inside the radial dis- 
tance r. (This formula is also used in the derivation of the 
Tolman-Oppenheimer-Volkoff equation.) In geometrized 
units we get for r = ro, 



4>'{r ) 



2m(r )/r§ + 87rr p(r ) 
l-2m(r )/r 



(20) 
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Since we are now concerned with what is going to be- 
come the throat at r = r , we need to evaluate po and 
P(q : o) m Eqs (fT2j) and (fl~3|) to obtain 

p(r) = p ( ro ) e ^(n,)(l+m)/2m e -<«r)(l+m)/2m ( 21 ) 



and 

p q (r) = e 2 ^[ Pq (r a ) - B]e- 2 ^ + B. (22) 
So Eq. (fTHI) becomes 



b'(r ) = 1 + 



167rr [-p(r )(l + 3m) - 2 Pq (r Q ) + AB] 
[2m{r )/r 2 + 87rr p(r )]/[l - 2m(r )/V ] 



(23) 

We already know that p — mp and that p q > 4B, so that 
the numerator is indeed negative. For the purpose of 
comparison, we will simply assume for now that Pq(ro) = 
AB. Then 



b'(r ) = 1 + 



167r[-/j(r )(l+3m) - 4B] 
[2m(r )/r 3 + 87rp(r )]/[l - 2m(r )/r ] ' 



(24) 

Disregarding r for the time being, suppose we determine 
m(ro) so that 2m(ro) exceeds 167r(4i3). We obtain 



16tt(4B) = 16tt(4)(145) 



MeV G 
(fin) 3 c* 



which leads to m(ro) ~ 1CP 8 m. The other term in the 
denominator, 87rp(ro), could be close to zero (since < 
m < 1). However, by simply doubling the value of m(ro), 
we can compensate for — p(ro)(l + 3ra) in the numerator, 
since this term is smaller than the other term. Moreover, 
taking p q to be 2 x 10 18 kg/m 3 , we see that in geometrized 
units, p q {ro) is in fact larger than and reasonably close 
to 4B, thereby justifying the assumption p g (r ) « 4B. 

Next, given the enormous density near the center of 
a neutron star, m(ro) 10~ 8 m would correspond to a 
fairly small value for tq: using p q (ro) ~ 4B again, 

^irrlp q {r ) = lO" 8 

yields ro«lm, although larger values would be allowed. 
So the above estimate holds. 

Finally, observe that the other factor in the denom- 
inator, 1 — 2m(ro)/ro, is now seen to be very close to 
unity. 

We conclude that b'(ro) < 1. Moreover, using Eq. (|22|) 
in Eq. (|17l) . we obtain b(r) — r < for r > r , so that 
6(r)/r < 1. 



to a thin-shell wormhole from a Schwarzschild black hole 
fl^il : while not part of the manifold, the black hole gen- 
erates the gravitational held. This comparison explains 
why our wormhole requires a minimum throat radius to 
help produce a sufficiently strong gravitational field. 



A. The need for quark matter 

If quark matter is eliminated from the model, then the 
right-hand side of Eq. (|14l) becomes 87rp(l + 3m). As a 
consequence, the 4S-term in Eq. ([24)) is also eliminated. 
Since 

|167r[-p(r )(l + 3m)] | > 2 x -vrp(r ) g^pfo), 

6 r 

the flare-out condition b'(ro) < 1 is no longer satisfied. 
It now becomes apparent that our conclusion depends 
critically on the fact that B is a constant: increase the 
throat size until 2m(ro) surpasses 16ir(4B). So based on 
our model, a topology change is possible only if there is 
a sufficient amount of quark matter in the center. More 
precisely, the assumption p(ro) « AB implies that the 
quark matter has to extend at least to r = tq. 



B. The interacting case 

We conclude with some brief comments on the interact- 
ing case. If the two fluids are assumed to interact, then 
the conservation equations take on the following forms 

(HOI: 



dp 
dr 



1 



2m 



and 



^ + 2^-5) 



p = Q 



-3Q. 



(25) 



(26) 



The quantity Q expresses the interaction between the two 
types of matter and falls off rapidly as r — ¥ oo. Further- 
more, since the interaction is assumed to be relatively 
weak, Q is a very small quantity compared to B. So in 
comparing m(ro) and p q {ro) hi the expression for b'(ro), 
the presence of the quantity Q is going to have little ef- 
fect. 



V. COMPLETING THE WORMHOLE 
STRUCTURE 

As we have seen, m(r ) is the mass inside the radial dis- 
tance r — tq. Since r = ro is now the throat of the 
wormhole, the interior r < r is outside the manifold. 
Although not part of the wormhole spacetime, it still con- 
tributes to the gravitational field. This can be compared 



VI. CONCLUSION 

Quark matter is believed to exist in the center of neutron 
stars. The analysis in this paper is therefore based on a 
two-fluid model comprising ordinary and quark matter 
with an isotropic matter distribution. It is shown that 
the extreme conditions may result in a topology change, 
that is, a neutron star could give rise to a wormhole. 
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